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Abstract 



We examine tiie possible tests of violation of the gravitational equivalence principle (VEP) 
at a muon storage ring via neutrino oscillation experiments. If the gravitational interactions of 



(N 

> . 

, the neutrinos are not diagonal in the flavour basis and the gravitational interaction eigenstates 

(N 



have different couplings to the gravitational field, this leads to the neutrino oscillation. If one 
starts with ^i^ beam then appearance of r^, e+ and /i^ in the final state are the signals for 
, neutrino oscillation. We have estimated the number of /i~ events in this scenario in — N 

\ deep inelastic scattering. Final state lepton energy distribution can be used to distinguish the 

\ VEP scenario from the others. A large area of VEP parameter space can be explored at a 

I ' future muon storage ring facility with moderate beam energy. 

: 

1 General Strategy 

X' 

■ The oscillation among the different neutrino flavours is now a well accepted solution to the solar 
neutrino problem [Q, ||. Recent results from the SuperKamiokande (SK) experiment @ at Japan 
on atmospheric neutrinos also support this proposition. The most popular explanation behind 
the neutrino oscillation is that neutrinos have non-degenerate masses and the mass eigenstates are 
not the same as gauge interaction eigenstates. Oscillation probability from one flavour to another 
is proportional to Am^, where Am^ is the difference between the square of the physical masses. 
Lots of work have already been done in this direction 

An alternative to this "mass-mixing (MM) solution" was proposed long back by Gasperini 
and independently by Halprin and Leung |5| . Their idea was that gravitational interaction of the 
neutrinos may not be diagonal in the flavour basis. Thus the gravitational interaction eigenstates 
are different from the weak interaction eigenstates. Now if these gravitational interaction eigen- 
states couple to the gravitational potential with different strength, neutrino oscillation can take 
place. This mechanism does not require neutrinos to have non-zero masses. 

The principle of equivalence has been the cornerstone of the general theory of relativity. 
Validity of this principle has been tested to a high precision for macroscopic bodies. But this has 
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not been tested experimentally in the microscopic and quantum regime. So it is prudent to keep 
an open mind towards this issue and try to prove/disprove the validity of equivalence principle in 
the current and future experiments. 

For a two flavour oscillation picture]^, the most general expression for the transition probability 
is given by, 

Pi^j = sin2 29 sin^ (^^^ (1) 

L, is the base-line length (distance which neutrinos traverse from the source to the detector). 6 
is the usual mixing angle. 

Expressions for A for the two cases (MM and VEP) are, 

A = ^ (MM) 

(VEP) (2) 



In the above equations, Ei, is the neutrino energy. A/ (= fi — f2) represents the difference 
between the coupling strengths of the two gravitational eigenstates with gravitational potential (j) 
and it quantifies the violation of equivalence principle in the neutrino sector. The basic difference 
between the two cases, on the dependence of E^, is evident. So their predictions can be different 
altogether. This alternative has attracted many people and a lot of work has been done |6|. 
Oscillation probability induced by VEP is, 

P = sin^ 29 sin^ (^2.538 x 10^^ LEyl^fcj^ (3) 

Here, L is expressed in kilometres and Ey in GeV . We will assume that (j) remains unchanged 

over the neutrino path. But, this is not very crucial for our purpose, we will parametrise VEP by 

the product (j) {= AF). Some recent work 0, analysing the SK solar neutrino data, gives the 

best fit value for this parameter. This comes out to be as small as 10~^^ for maximal mixing. 

In this letter, we will try to explore the possible signals of VEP at a muon storage ring. Here, 

we assume the sample design, for muon production, capture, cooling, acceleration and storage 

as given in ref. |Q. Number of available muons directed to the neutrino detector of per year is 

10^0 (one can have //^ or ) . For the purpose of illustration, we start with beam. These 

positively charged muons will produce ~ 10^*^ z^^s and same number of fgS. Muon anti-neutrinos 

(or electron neutrinos) thus produced, traverse a distance L before colliding on a fixed target. For 

a sufficiently energetic initial muon beam, neutrino-nucleon (in the target material) interaction 

is in the deeply-inelastic regime. Now, on their way to the detector, if some P^s are oscillated to 

v'e or Pt-s, in the final state we may observe rs or es with the same sign with the initial muon 

beam, via neutrino nucleon charged current interaction. The fgS (also coming from fj, decay) may 

transform to f^s or z^t-s, which can be the source of fi~s or r~s. We will call them wrong sign 

fis or rs, as they have opposite charge to the initial ^-beam. The un-oscillated (z^e)s can only 

^we will stick to this in the present analysis due to simplicity. But this is sufficient to illuminate the underlying 
principle. 
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give rise to /i"*" (e~)s. So appearance of r leptons of either sign and negatively charged muons 
or positively charged electrons are definitely the signal for neutrino oscillation at a muon storage 
ring with positively charged muon-beam. These issues have been studied in detail earlier |^] in 
the context of MM scenario. Here we will mainly concentrate on the results for fi~ appearance in 
the final state. This means our focus is on the i^e ^ '^fi oscillation due to VEP. We will see in the 
following, that fi~ event rate in MM scenario is negligibly small with presently SK allowed values 
for Am-g^ and sin^ 26. So appearance of considerable number of fi^ events at a muon storage 
ring, is a signature which goes against the MM solution to neutrino oscillation problem. So the 
charge identification for the muons is necessary for the study of neutrino oscillation physics. This 
seemed to be achievable according to the ref . |^ . VEP solution does not fit very well to the — Ur 
oscillation data from SK. So we will not calculate the r+ appearence rate in the following. 

For lower values of energy of the initial muon beam, r~ cross-section is slightly lower than the 
fi~ case due to the phase space suppression. Furthermore, for tau in the final state, above cross- 
section has to be multiplied by appropriate branching ratio of r. Generally one prong hadronic 
or leptonic decay channel is the cleanest way to detect a r. Former has a branching ratio of 45 % 
whereas the later has of 17 %. 

Number of fi~ events coming from neutrino-nucleon DIS can be obtained by folding the charged 
current cross-section by oscillation probability, neutrino flux and finally by the total number of 
nucleons, A'^, present in the fixed target. This last quantity depends on the amount of the target 
material present. 

N,, = NnJ aiu^ + N^fi-+X)^ P{u, ^ u^) dE,^ (4) 

We present our results for 1 kT of the target material^. For any other detector size these results 
can easily be scaled. The charged-current neutrino-nucleon cross-section is calculated assuming 
almost the same number of protons and neutrons, present in the target material. The expression 
for the same is not given here but can be obtained in several places [^, ^. Our results are based 
on a simple parton level monte-carlo event generator. We have not incorporated any detector 
effects. Even incorporating these, the essence of our analysis will remain the same. We have 
used CTEQ4LQ parametrisations |jl^ for the parton distribution functions to estimate the above 
cross-section. 

Neutrino fiux depends on the number of neutrinos, initially present. At a muon storage ring 
this number is equal to the number of muons present in the beam. Neutrino flux also depends 
on energy/angular distribution of the neutrinos coming from ^ decay. We will elaborate more on 
this in the following. 

The area of the target and the base-line length, L, define a cone with half angle Od, which can 
be written as Rd = L Od- {Rd is the radius of the detector if we assume it has a circular cross- 
section.) Thus the choice of detector size for a long-baseline experiment must differ with that of a 
short base-line experiment and also depends on the energy of the muon beam. If the muon-beam 
^This corresponds to N„ = 6.023 x 10^^ 
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energy is increased, neutrinos coming from muon decay, become more and more boosted in the 
muon direction. This effectively increases the luminosity of neutrino beam. A higher effective 
luminosity of the incident neutrino beam can be achieved with a sufficiently high energy beam 
and /or by decreasing the effective area of the target (keeping the amount of target material fixed) . 
In our event generator, we only consider those us, for which O^^i, < O^- is the angle between 
the decaying muon beam and the neutrino.) Of course this is fully determined by the muon decay 
kinematics. In this analysis we will assume a fixed detector size with R4 ~ 6 mts.. Target area 
corresponding to this radius are similar to that of proposed ICANOE detector 



2 Discussion of the Results 



In Fig. H, we present variation of the number of wrong sign muon events coming from i^^j-N DIS 
experiments for VEP scenario. Number of tau events, from v^--^ scattering will almost be the 
same as this. We have chosen AF = 2xl0~^^ and sin^ 26 as 1. Choice of these parameter values 
are from ref. Q. They have used the SK solar neutrino data, to calculate the best fit values of 
/S.F and sin^ 29. The only other bounds on these two parameters come from the LSND and E776 
experiment at BNL |12]. Ref. [0] constrains AF responsible for the *^ t^x {x = fi,T or stetrile) 



oscillation. The other one [12| constrains the same parameters for ^ i/^ oscillation. 
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Figure 1: Number of ^ events coming from u^-N DIS at a muon storage ring with AF =2 xlO 
and sin^ 29 = 1. Ad is the area of the target. 



We will not present here the event rate from Vr—N scattering. VEP solution to atmospheric 



Vfj, — i't oscillation is disfavoured by the the SK data [0, In ref. |15, |T^, a analysis 



has been done with SK atmospheric neutrino data, assuming a general power law dependence of 
oscillation probability on the neutrino energy (~ E^). Minimum is obtained for n = —1 i.e. 
for the MM scenario. 

But it should be borne in mind that along with this indirect evidence (x^ analysis in the ref. 
|T5| , ) against the VEP solution to the atmospheric — oscillation, one also should look for 
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some direct evidences to disprove such propositions. Neutrino factory with a huge and precisely 
known neutrino flux is an ideal place for this, z/^ <-> Uj- oscillation due to VEP and subsequent 
Ur — N scattering can produce r"*" events. The search for r appearence in final state will be one 
of the major goals of neutrino oscillation experiments. So we would not have to pay some extra 
price to look for signal which can can arise from VEP in neutrino sector. 

We have checked that event rates remain almost unchanged for two baseline lengths (First 
one, 250 Km., corresponds to a baseline from Kamioka to KEK and the second one, 740 Km., 
corresponds to a baseline from CERN to Gran Sasso or from Fermilab to Soudan experiment.) 
over the energy range we used in fig. ^. For the values of neutrino energy and AF we have 
used, oscillation probability grows quadratically with L. At the same time, for a fixed target area 
and size, effective luminosity of colliding neutrinos decreases with L, due to the decrement of the 
detection cone subtended by the target. These two effects compensate each other, for the detector 
area we have chosen. We checked that for a sufficiently large target area, and/or sufficiently high 
muon beam energy, when almost all the neutrinos coming from muon decay can be intercepted 
by the target, (i.e. when neutrino flux becomes independent to L) the event rate indeed grows 
quadratically with L. 

Number of /i~ events in MM scenario, when calculated using SK allowed values for mass^ 
difference and mixing angle (Am^^ = 2 x 10~^el^^, siv?29 = 1. These values are for MSW 
solution. Vacuum solution gives Amg^ ~ lO^^'^eF^ [0)' comes out to be small compared to the 
VEP scenario over the range of muon energies we are interested and is beyond the scale used in 
fig.||. The other source of /u~ is from the r~ (produced in z/,- — interaction) decay. But this 
rate is further suppressed by r ^ branching ratio. Thus the appearance of ^~ in final state 
(starting with fi^ beam), is a definite signature of neutrino oscillation not originated by MM 
scenario. Appearance of /U~ in the final state may also be accounted by the R-parity violating 
/x-decay and/or v — N interaction [19|]. 

We confine our discussion to the appearance of wrong sign muons in the final state. Wrong 
sign r event rate will be somewhat lower (due to phase-space suppression) than this if difference 
of couplings between i/^ and v^- to the gravitational potential is at the same ballpark of i/^ and Ve- 
Apart from ref. Q, the only bound on the AF, coming from two terrestrial experiments (LSND 
and E776) are not consistent with each other. And these values are order of magnitude greater 
than the best fit value obtained in ref. Q. Upper bound coming from the E776 experiment is 
around 3x10^^^ for sin^ 2^ = 1. On the other hand LSND results predicts a lower bound on 
AF (~ 5 X 10^^^) for the same value of the mixing angle. The SK situation is a bit different 
from the other two cases (In terrestrial experiments, neutrinos are mainly affected by the earth's 
gravitational field. In the former, neutrinos traverse the sun-earth distance. Gravitational fields 
of sun and near by stars along with the earth have to be considered in such a case.). We want 
to stress that even with such a small value of the above parameter, number of events in VEP 
scenario will be much larger than that of MM scenario. In the following we will explore the 



^R-parity violating theories allow lepton flavour violation, thus /i or r appearance in such a theory can also be 
explained without flavour oscillation. 
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possible region in AF — sin^ 26 space for which a perceptible signal can be obtained such that we 
can exclude that region at a certain confidence level. 
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Figure 2: Normalised energy distribution at a 50 GeV muon storage ring. Solid curve is for 
MM scenario and the dashed curve is for VEP scenario is the target area. 

Before we go into the exclusion contours, we want to point out another distinguishing feature 
of the VEP scenario, which is also an artefact of the energy dependence of oscillation probability. 
We look at the normalised energy distribution of the /i~ in the final state. In MM scenario lepton 
energy distribution falls of rapidly in contrast to the VEP scenario. This is reflected in fig. |2|. 
If the VEP parameter /S.F is so small that event rate is comparable with event rate in the MM 
scenario ^, energy distribution of the final state lepton can be used to discriminate between these 
two scenarios. So one can have an idea about the energy dependence of the oscillation probability 
from the lepton energy distribution. If one concentrates on the /i~ appearance rate, then MM 
scenario is outnumbered by VEP scenario. Strong constraint can be put on this scenario from 
the wrong sign muon search. We have already mentioned about /x" or r+ appearance due to 
the R-parity violating interaction at a muon storage ring. Wrong sign muon or tau events can 
be accounted without the oscillation phenomena in such a kind of theories. Here also the lepton 
energy distribution is expected to be different from MM scenario. 

Until now, we have used the best fit values for IS.F and sin^ 26 as given in ref. ||^], but the 
main purpose of this letter is to point out that even a tiny violation of equivalence principle can 
be detected at a muon storage ring facility and can be distinguished from the MM scenario quite 
efficiently. In the following we like to show, in the AF — sin^ 26 plane, the region which can be 
excluded at 95 % C.L. 

Prom our previous discussions, it is clear that /i~ appearance search is the best probe for 

VEP at a neutrino factory. Firstly, //~ event rate is higher than the r~ event rate. And with 

*One can see from previous discussion, /x~ event rate in MM scenario is negligibly small compared to the same 
in VEP scenario with the values of parameters we have used. But having larger amount of target material and at 
the same time smaller value of AF can make these two comparable with each other. 
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experimentally allowed values of MM parameters, one cannot have /i~ signal from MM scenario. 
We will present the exclusion contours from fi~ appearance search for two values of beam energy. 
We will not discuss here in details the possible backgrounds for such signals. These issues have 



been discussed in ref. |2C]. We have applied the same kind of cuts to reduce the background. 



Pt 



> 2 GeV, 



AR^-x > 0.4 



(5) 



AR^-x is the isolation between the ^~ and the DIS hadronic products. After applying the above 
set of cuts the signal efficiency comes out to be 70 %. The other source of wrong sign muon is 
from the decay of the wrong sign taus' coming from the Vr — N scattering. These come from 
fg <-> oscillation. So this is also due to flavour oscillation and add to our signal. But we confine 
ourselves to two flavour oscillation only and do not add the numbers of wrong sign muons coming 
from tau decay. 
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Figure 3: 95 % exclusion contour in AFg^ — sin^ 29 plane for 50 and 20 Ge V muon storage ring 
from wrong sign muon signal. 

The exclusion contour depicts the fact that for small values of AF one need maximal mixing 
between the flavours. For higher values of VEP parameter very small mixing would sufficient to 
produce enough number of signal events. For AF less than 10^^^ , sin^ 29 and AF are correlated 
in a linear fashion. In this region, oscillation wave-length A (see eqn. |l]) is much larger than 
the base-line length, L. For larger values of AF , L is comparable with A. Here, siv?{AF Ey L) 
attains its maximum value. This explains the correlation between the two oscillation parameters 
near the 'knees' of the contours. For even higher values of AF, L >> X and the sin^(^) can 
be approximated by ^. This causes the sharp rise of the contour almost independent of mixing 
angle. 

A different phenomenon which also cause the neutrino oscillation with the same neutrino 
energy dependence of oscillation probability, is the violation of special relativity (VSR) [^] . One 
can parametrised the oscillation probability in the same way. So the above analysis is completely 
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applicable to VSR also. The energy distribution and exclusion contours derived above would be 
the same for the later. 

3 Conclusions 

To summarise, we have shown that at a muon storage ring with moderate energy, neutrino flavour 
oscillation due to VEP can be probed. The appearance event rate due to u^ — N DIS has been 
calculated within RG improved parton model. We compared the VEP results with the MM results 
and have shown that, the former will give overwhelmingly large number of wrong sign muon events 
than the later. In VEP model, the number of /i~ events increases with beam energy. The energy 
distributions of the final state muon have also been compared for two different scenario. If the 
violation of equivalence principle is so small that it would produce number of events comparable 
with that of mass mixing model, then energy distribution of the final state lepton can be used as 
a discriminator. Finally we have shown the 95 % C.L. exclusion contours in the AFe/^ — sin^ 26 
plane for two different muon beam energy. 
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